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We derive the complete form of the van der Waals dispersion interaction between two infinitely 
long anisotropic semiconducting/insulating thin cylinders at all separations. The derivation is based 
on the general theory of dispersion interactions between anisotropic media as formulated in [J. N. 
Munday, D. Iannuzzi, Yu. S. Barash and F. Capasso, Phys. Rev. A 71, 042102 (2005)]. This for- 
mulation is then used to calculate the dispersion interactions between a pair of single walled carbon 
nanotubes at all separations and all angles. Non-retarded and retarded forms of the interactions 
are developed separately. The possibility of repulsive dispersion interactions and non-monotonic 
dispersion interactions is discussed within the framework of the new formulation. 

PACS numbers: 78.20.Bh, 34.30.-h, 77.22.-d 
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I. INTRODUCTION 



Single walled carbon nanotubes (SWCNTs) constitute a unique class of materials with chirality-dependent dielectric 
properties^ that give rise to interesting consequences in terms of their van der Waals - dispersion interactions. Several 
experimental procedures have exploited the differences among these properties in order to separate SWCNTs by 
chirality (see Ref. |2| and references therein). In order to reliably separate a polydisperse solution of SWCNTs into 
monodisperse fractions of one chirality, one needs to understand the detailed features of the interactions between 
constituent SWCNTs. Apart from this, the dispersion interactions could also drive micro and nano-mechanical 
actuators that could transduce rectilinear motion or even convert rectilinear into rotational motion 3 . To achieve 
• this goal a much more rigorous understanding of fundamental forces and, in particular, van der Waals - dispersion 
| interactions between SWCNTs is needed. Several important advances have already been made^. Our work continues 
. these efforts. 

We derive the complete van der Waals - dispersion interaction free energy between two anisotropic cylinders at all 
separations, including the retardation effects. We start with the Lifshitz theory of van der Waals interactions between 
two semiinifinte anisotropic uniaxial dielectric layers across a finite layer of dielectric function e m and thickness I 
as worked out by Barash 5 , the result of this calculation is the interaction free energy between the two layers as a 
function of their separation I and the angle between their principal dielectric anisotropy axes 9: Q(£, 9). The dielectric 
response of the two dielectrically uniaxial half-spaces is given by the values of their dielectric functions ej[, parallel 
and e±, perpendicular to their respective axes. We shall use eTJ (eTjZ) and eTJ[ {WX) f° r the left and right half- 
spaces, respectively. Note also that in the theory of van der Waals interactions 7 - all the dielectric response functions 
are evaluated at imaginary frequencies, thus en j_ = en j_(iu>). en is referred to as the London - van der Waals 

transform of the response function en and is given by the Kramers - Kronig relations. It is strictly a real, 

monotonically decaying function of to. 

From the interaction free energy between two half-spaces one can extract the interaction between two cylinders by 
assuming that the two half-spaces arc dilute assemblies of anisotropic cylinders. This derivation closely follows the 
arguments of Pitaevskii for evaluating the interactions between isotropic impurity atoms in a homogeneous fluid 8 . We 
assume that the two anisotropic half-spaces are composed of anisotropic cylinders at volume fractions V\ and i>2, with 
£ c i,_l (c c 2,_l) and e c 1 ^ (e C 2,||) as the transverse and longitudinal dielectric response functions of the cylinder materials. 
We then expand Q{i,ff) for two half-spaces as a series in v\ and V2 and evaluate the coefficient multiplying the V1V2 
term. The volume fractions v\ and t>2 scale with the area density of the cylinders (N\, N2) in the direction of their 
long axes as v\ — N\ ttR\ (v2 — N2 ttR^)- It then follows^ that the interaction free energy between two cylinders, 
G(£,9), whose axes are contained within the two parallel boundaries at a separation £, but skewed at an angle 9 is 
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given by 

d 2 G(£,9) 
d£ 2 



NxN 2 sin 9 G(£,6). (1) 



Conversely, the interaction free energy per unit length, g(£), between two parallel cylinders is given by the Abel 
transform 



= N,N 2 / gW£ 2 +y 2 ) dy. (2) 



d 2 g(e,9 = o) 
dp 

In both cases we expand Q(£, 9) to find the coefficient next to viv 2 (or equivalently NiN 2 ), take the second derivative 
with respect to £, then use Eqs. [1] and [2] in order to obtain the appropriate pair interaction free energy between 
cylinders. Note that such an expansion is possible only if the dielectric response at all frequencies is bounded. In 
the case of an ideal metal Drude-like dielectric response this expansion is not feasible and our method can not be 
transplanted to that case automatically. 

The closest attempt in the literature to evaluate the interaction between two cylinders at all separations comes from 
Barash and Kyasov 9 . Where this approach can be compared with the one presented here, i.e. for two parallel isotropic 
cylinders, the results for the interaction free energy between parallel cylinders coincide completely. As far as we are 
aware our calculation thus presents the first attempt to evaluate the van der Waals - dispersion interaction between 
two skewed anisotropic cylinders at all angles and at all separations. Its major drawback is that by construction it 
is valid only for materials with finite dielectric response and thus cylinders with ideally metallic dielectric functions 
showing a Drude-like peak at zero frequency can not be treated by the theory presented below. Also the calculation 
is only valid for infinitely long cylinders and finite length effect are not taken into account exactly. 

II. DERIVATION 

We use the Pitaevskii ansatz in order to extract the interactions between two infinite anisotropic cylinders at 
all separations and angles from the interaction between two semi-infinite half-spaces of anisotropic uniaxial dielectric 
material. We start with the fully retarded van der Waals - dispersion interactions between two semiinfinite anisotropic 
dielectric slabs^. The full interaction form is quite involved, but it has a simple limit if the two semiinfinite slabs, 1Z 
and C, separated by an isotropic medium of thickness £, are composed of rarefied material. 

In order to get the interaction free energy between two anisotropic cylinders we assume that both semi-infinite 
substrates (half-spaces), C (1) and TZ (2), are composite materials made of oriented anisotropic cylinders at volume 
fractions v\ and v 2 , with e c i,j_ (e c 2 ,j_) and e c x . 1 1 ( £C 2,||) as the transverse and longitudinal dielectric response functions of 
the cylinder materials. For the semi-infinite composite medium of oriented anisotropic cylinders with local hexagonal 
packing symmetry, so that the corresponding cylinder volume fraction is v, the anisotropic bulk dielectric response 
function can be derived in the form (see Ref. @, p. 318) 

ej[ = e m (l + i;A||), e± = e m f 1 + TZ^j ) ' ( 3 ) 
where the relative anisotropy measures in the parallel and perpendicular direction are given by 

A ± = t±_ ±1 A|| = tlZ^L. (4 ) 

C _L ~r Cm Cm 

In our case, this holds for both C and TZ half-spaces with the appropriate volume fractions and dielectric responses. 
e m is the dielectric function of the isotropic medium between the cylinders as well as between regions C and 1Z. We 
assume in what follows that all the response functions are bounded and finite. 

The formulae in Eqs. 1 1121 connect the interaction free energy of two semiinifinitc half spaces with the interaction 
free energy between two cylinders either parallel or skewed at a finite angle 9. The Barash result^ for the complete 
retarded form of the interactions between two uniaxial media, G(£,9), is quite complicated (note also a typo that 
propagated starting from the original version of the calculation^) but can be straightforwardly expanded to second 
order in N (a term proportional to V\v 2 ) for the dielectric response functions of the form Eq. [3l yielding the following 
result 

d 2 g(£,9) k B r^' r n , n ffm 

Q d Q Jn ■ ( 5 ) 



d£ 2 2vr ^ n ^ ^ dP 

n— 
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In the above equation, n represent the (thermal) Matsubara indices, the prime on the summation means that the 
weight of the n — term is 1/2 (see Refs. l7||9j for details). The second derivative of the function f(£,9) can be 
obtained explicitly in terms of the ratios between the relative anisotropy measures (Eq. |4|) defined as 

a = —a = 2~. T7 r (6) 

A|| (e c ± + e m )(e c || - e m ) 

and is obviously frequency dependent. Parameters a\ and a 2 can be thought of as a specific measure of the anisotropy 
of the cylinders in the left and right half-spaces when compared with the isotropic bathing medium m. Note that 
they vanish when the transverse dielectric response of the cylinder material equals the medium response. The explicit 
form of the second derivative of f(£, 8) now follows as 



d 2 f(t,9) 



V\V2&i 11 A 2 II e 



32 



(Q 2 



,-2 ) 



(1 + 3ai)(l + 3a 2 )Q 4 + 2(1 + 2a x + 2a 2 + 3aia 2 )Q 2 e m ^f + 2(1 + ai)(l + a 2 )e m 2 ^ 



+ (l-ai)(l-a 2 ) (Q 2 + 2e 



"J 7: 



cos 26 



+ 



(7) 



Here R\ and R 2 are the cylinder radii, assumed to be the smallest lengths in the problem^. The frequency dependence 
of the dielectric functions is in e m {iuj n ), e c ±(iuj n ) and e c y («tt>„), and therefore also a — a(iuj n ). The frequencies in the 
Matsubara summation are co n — 2tt ^jf-n. Note that Eq. [7] is symmetric with respect to 1 and 2 indices (left and 
right half-spaces), as it should be. 

This is as far as a general theory can go. We must now deal separately with the cases of skewed and parallel 



d 2 g(e,9) 
IP 

two cases, see Eqs. 11121 We first analyze the case of skewed cylinders. 



cylinders, since the connection between a y d p' a> and the effective pair interaction between cylinders is different for the 



A. Skewed cylinders 



We use Eq. [T]to obtain the interaction free energy between two skewed cylinders: 

G *' (9 > = ~T~A Ti ■ a Z-v A iJI A 2,ll / udu - ^ 2 s g(a 1 ,a 2 ,u,p n ,6), 

64-7T rsmff ^— < J n (ir + pi) 

n— 

where u = Q^, 



g(ai,a 2 ,u,p n ,6) = 2 [(1 + 3ai)(l + 3a 2 )u 4 + 2(1 + 2a x + 2a 2 + 3aia 2 )u 2 p 2 + 2(1 + ai)(l + a 2 )j4] + 

+ (l-a 1 )(l-a 2 )(u 2 + 2p 2 n ) 2 cos29 (9) 

2 

and p 2 = e m («w r j)^f-£ 2 . Another change of variables with u — p n t, yields 

l. ti 2p2p2 00 ' roc -2p„Vi 2 + l 

G(^0)--^| r7I ^fE A i.ll A 2.llPn / ^ (t2 + 1) g(t,oi(w w ),oa(ia; TO ),g), (10) 

n— ^ 

with 



<7(i, a 1 ,a 2 ,9) = 2 [(1 + 3ai)(l + 3a 2 )i 4 + 2(1 + 2a x + 2a 2 + 3 ai a 2 )t 2 + 2(1 + ai)(l + a 2 )] + 

+ (l-a 1 )(l-a 2 )(t 2 + 2) 2 cos26». (11) 
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This is the final result for the cylinder-cylinder interaction at all angles when the radii of the cylinders are the 
smallest lengths in the system. It includes retardation and the full angular dependence. Some simple limits that can 
be obtained form this general expression. 

The non-retarded limit where c — ► oo, has already been explored in Ref. 0. There p n — ► for all n and we obtain 
from Eq. 1 



G(£,9;c — > oo) 



^Ai,||A 2 ,|| / u 3 du e~ 2u [2(1 + 3a 1 ){l + 3a 2 ) + {l-a 1 )(l-a 2 ) cos 26] 



k B T n 2 R 2 R 2 
~ 64tt £ a sin 9 

U rp 2 r>2 r>2 00 ' o 

= —^l^dT, A ^ll A 2-ll g + 3a ^ + 3a2 ) + 0- - - «2) cos 29] . (12) 

n=0 

This formula could also be obtained directly from Eq. [10] taking into account that in the t integration only the terms 
with large t contribute to the final integral. Expanding the whole integrand for large t returns us to Eq. Q21 

At low temperatures, when the summation over the Matsubara frequencies can be turned into an integral over n 
with dn — h/(2n kBT)dw, the corresponding interaction free energy is 



* n 2jl2 R 2 roc foo e ~-2y/e m (iuj)feVf 2 + l 

G{L ' 6) = ~ 128tt2 dhko J dW w4A i.llMA 2 ,|| Me m M 2 ^ tdt g{t, ai {™)M™)>0)- ( 13 ) 



We now rework this equation to obtain the retarded result for the interaction between two semiconducting cylinders. 
Note here that we can not derive the Casimir limit properly as our formulation is not valid for nominally infinite 
zero-frequency (Drude-like) dielectric response. For that case see Ref. & First instead of variable lo, we introduce 
x = -\/t 2 + 1 uj. Then, following closely the arguments in Ref. [Ill we obtain the interaction free energy in the form 



GM= - 12W kmf £t " (Q) Al '" (Q)A2 '" (0) I dxX J Q ^ 1 f2—y^9(t,a 1 (0),a 2 (0),9). (14) 

Here e m (0) and ai(0),a 2 (0) denote the static, i.e. zero frequency, values of the corresponding functions. Obviously 
in this regime the interaction free energy decays faster with separation, being a reflection of the retardation. All the 
frequency dependence of the material properties is reduced to the static response in this limit, just as in the general 
Lifshitz analysis^. 

B. Parallel cylinders 

The analysis here is somewhat more complicated because the pair interaction energy between the cylinders involves 
the inverse Abel transform. We start with 

d 2 g(£,9 = 0) k B T^' f°° „.„<Pf(i,0 = 0) 



d£ 2 2ir 4- ' J ^ * d£ 2 



71=0 



where 



d 2 f(£,9 = 0) _ wi«aAi f „A 3 ,|| e 2, \ /( ' J+< " 



d£ 2 32 



(Q 2 



,2 4 " 

(1 + 3t»i)(l + 3a 2 )Q 4 + 2(1 + 2oi + 2a 2 + 3 ai a 2 )Q' A e m ^ + 2(1 + oi)(l + a 2 )e m 2 ^ 



+ 



+ (l- ai )(l-a 2 )(g 2 + 2e m ^)^. (Hi) 



and again V\ = N nR 2 (« 2 = N 7ri? 2 ) and a — We continue by introducing the Abel transform and its properties. 

Namely, if we define 



+oo 



gWP + V 2 ) dy = f(y), (17) 
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then 



9il) =-ir miL. ,18, 

TT J, jy2 _ P 



Taking this into account when considering Eqs. 1161 we remain with 



sW = — %T R i R V>2 A i,|| A 2,|| / / „ a0 \ Q d Q— : —2—h(ai(^n),a 2 (iuJn),Q,£ m (iuJn 

61 n Jl \/W ! -i z Jo (Q 2 + e m (iu„ )-*Vi 2 



where 

h(a 1 ,a 2 ,Q,e m -^-) = 2 



2 " 

(19) 



2 4 

(1 + 3o x )(l + 3a 2 )g 4 + 2(1 + 2ai + 2a 2 + 3a 1 a 2 )Q 2 e m '^ + 2(1 + oi)(l + a 2 )e m 2 ^ 



(l- ai )(l- fl2 )(g 2 + 2e m 4) 2 - (20) 

cr 

2 

As before, wg introduce p n — c m (iw n )^f^ 2 , m = Q£ and ?/ — ► yjl. This allows us to rewrite the above integrals as 
k B T R\R\ . . /"+ 00 Z" 00 , e - 2!, V u2+p » v ,. v 2^ 

9 ( > = Al .H A ^ll _/ "TTTl y o ( u 2 +p 2)l/2 M°l(* w n),02(lWn),« > P„), (21) 



n=0 

and 



h(ai(iuj n ),a 2 (iuj n ), u,p 2 n ) = 2 [(1 + 3oi)(l + 3a 2 )u 4 + 2(1 + 2ai + 2a 2 + 3aia 2 )u 2 p 2 + 2(1 + oi)(l + a 2 )p£] + 

(l- ai )(l- a2 )( M 2 + 2p 2 ) 2 . (22) 

This is the final result for the interaction between two parallel thin cylinders at all separations and contains retardation 
effects explicitly. In general, the above expression can only be evaluated numerically once the dielectric spectra of 
component substances are known. 

In the non-retarded limit, c — > oo, the above formula reduces to 

U J 1 c>2 p2 °° ' r+oo 7 ,.00 

fl (4c— oo) - Ai.iiAa.nfa + Sfax + aaJ + igoiOa) / \ u*du e" 2 ^ = 

^ * rr; ji vr - 1 Jo 



^2 d2 D 2 . 



9 fc B T n R R 

- - — ^-^y^ A 1: y A 2 ,| {3 + 5[ai(iw„) + a 2 (iw„)] + 19ai(iw n )a 2 (iw„)} . (23) 

For the case where the two interacting cylinders are composed of solid isotropic dielectric materials this form of the 
interaction free energy can be compared with the result obtained by Barash and Kyasov ( Eq. 10 ir>^) and can be 
reduced to it exactly. 

As with skewed cylinders, we can take the zero temperature limit where the summation over the Matsubara 
frequencies becomes an integral over n with dn — h/(2n fcsT)dw. Again we introduce x = ^y/W+l to. Then, as for 
skewed cylinders, we obtain the interaction free energy per unit length of two parallel cylinders, 

*_ _2 p2 p2 roo r+oo j roo ,j, -2-i/e m (0) yx 

3( ^-6&^ £ ™ (0)5/2A ^» (0)A ^' (0) l dxx5 l V^=U (, 2 + 1)7 /2 MtM0)M0)). (24) 
Here 

h(t, oi, aa) = 2 [(1 + 3oi)(l + 3a 2 )i 4 + 2(1 + 2a x + 2a 2 + 3 ai a 2 )t 2 + 2(1 + oi)(l + oa)] + (l-ai)(l-a 2 )(i 2 +2) 2 . (25) 

The spatial dependence is, again, one power higher in the retarded regime than in the non-retarded regime. All the 
frequency dependence of the material properties in the retarded limit is again reduced to the static response as in the 
Lifshitz analysis^ 1 -. 
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III. NUMERICAL RESULTS: RETARDED VS. NONRETARDED INTERACTION 

Expressions for the van der Waals - dispersion interaction free energy between parallel and skewed cylinders derived 
above can be analysed numerically, once the spectral properties of the interacting cylinders are given. In all the 
results in this section, the medium between cylinders is assumed to be vacuum, 63(^0;) = 1, and the temperature is 
T = 297 K. The spectral properties of the cylinders are taken to be those of thin single walled nanotubes as calculated 
by ab initio methods in the optical ranged. Robust quantum mechanical codes have been developed to give us the 
very accurate spectral data needed to investigate the numerical consequences of the above theory. Without going 
into details, we should note that the dielectric spectral data depend on chirality and some SWCNTs even exhibit 
significant optical anisotropy between their radial and axial directions^. For purposes of illustration of the theory 
developed here, we choose [5,1] and [29,0] semiconducting carbon nanotubes, which have been previously analyzed 
and which differ substantially in their radii as well as spectral properties^. This choice is motivated by the fact that 
although [4,2] is technically the smallest SWCNT - it has a radius of only 0.207 nm - the [5,1] and [29,0] have a bigger 
difference in their spectra and are thus more appropriate to explore the effect of disparate dielectric spectra on van 
der Waals - dispersion forces. 

The theory developed here should be safely applied to this case since the tubes are not metallic and their dielectric 
response functions are thus always bounded. Furthermore, their dielectric responses and radii (1.135 and 0.218 nm 
for [29,0] and [5,1] SWCNTs, respectively, when measured from the center of the cylinder to the centers of the carbon 
atoms) are very different, which is another reason for choosing these tubes as the benchmark for application of the 
theory. For illustrative purposes we also disregard the finite core size of the [29,0] SWCNT that would require a more 
careful modeling of its effective dielectric response 2 - and thus introduce additional parameters that would complicate 
the understanding of the retardation effects in van der Waals - dispersion interactions between these two SWCNTs, 
which is our primary aim in this paper. Once the surface to surface separation between two SWCNTs is greater than 
approximately two SWCNT outer diameters 2 this approximation turns out to work quite well. 

A. Parallel cylinders 

It seems plausible that retardation effects are largest in the parallel configuration. Although the calculations are 
most demanding in this case, we nevertheless analyze it first. Figure [1] shows the interaction free energy per unit 
length for two parallel [5,1] SWCNTs (panel a), two parallel [29,0] SWCNTs (panel b) and parallel [5,1] and [29,0] 
SWCNTs (panel c). 

It is obvious, Fig. [TJ that the departure of the retarded results from the nonretarded ones takes place only when 
separations are larger than about 50 nm, but the strength of van der Waals interaction there is weak, on the order of 
10~ 8 ksT I nn*i£ for [5,l]-[29,0] pair of SWCNTs. One can see very gradual deviation of the retarded results from 
£~ 5 to £~ 6 scaling with the separation distance (note log-scale on both axes). 

For £ = 2 nm, the fully retarded value of the van der Waals interaction for the [5,l]-[29,0] pair is -0.13294 ksT 
I nm, while the nonretarded value is -0.13377 ksT / nm, so the contribution of retardation at this distance is only 
about 0.6 %. The nonretarded value of interaction then is of some use for a large interval of intercylinder separations 
due to the strict power-law behavior of the nonretarded values cx £~ 5 . For example, for £ — 16 nm, the interaction 
is -0.13377 (2/16) 5 = -4.0823 10~ 6 fc^T / nm. The full retarded expressions give a mostly negligible contribution to 
the van der Waals interactions at spacings of relevance. Its importance can be boosted when considering interactions 
between different types of SWCNTs in a dielectric medium with appropriate (or even "tailored") response. In that 
case and depending on the full dielectric spectra of the two cylinders and of the medium, retardation coupled to the 
dielectric spectra can in principle lead to a change in sign of the interaction via a similar mechanism first considered 
for the interaction of ice with vacuum across a liquid film, as studied by Elbaum and Schick^. (In that case, the 
two half-spaces are ice and vacuum, while the dielectric medium between is water). These effects will be investigated 
in Sec. IIVI Note that the sign change can also be obtained for hollow cylinders at small separations, but it has a 
different origin^. 

Of all the three cases considered, the van der Waals attraction is the strongest for the pair of [29,0] SWCNTs 
(-0.5325 fc^T/nm at £ = 4 nm). The most important reason for this is the explicit dependence of the van der Waals 
interaction on radii of the two cylinders (RfR%). The radius of [29,0] SWCNT is 5.21 times bigger than the radius of 
[5,1] SWCNT. 
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FIG. 1: The van der Waals interaction free energy per unit length between two parallel [5,1] SWCNTs (panel a), two parallel 
[29,0] SWCNTs (panel b), and parallel [5,1] and [29,0] SWCNTs (panel c) as a function of their separation. Circles represent 
the nonretarded expression^, while squares represent the fully retarded calculation Eq. 1211 Dashed and dotted lines indicate 
£~ 5 and £~ 6 dependences, respectively. The interaction free energy is calculated only for interaxial separations larger then 
I = R\ + i?2, i.e. in the non-intersecting regime of the cylinders. Close to this value the thin cylinder approximation breaks 
down and is superseded by a different type of calculation, see^ for details. 



B. Skewed cylinders 



These calculations are numerically less demanding since the dimensionality of integration is one rather then two 
for parallel cylinders. Figure presents the dependence of the van der Waals interaction free energy for two [5,1] 
SWCNTs (panel a), two [29,0] SWCNTs (panel b) and [5,1] and [29,0] SWCNTs (panel c) crossed at the right angle 
(8 = 7r/2). Again, the contribution of retardation is small. Its observable effects in this configuration take place at 
about 30 nm, similar as in the case of parallel nanotubes. The nonretarded and the retarded values for the van der 
Waals interaction for [5,l]-[29,0] pair separated by t = 2 nm are -0.3255 fcgT and -0.3285 fegT, respectively. The 
contribution of retardation effects to the van der Waals interaction at this distance is thus 0.9 %, similar to the case of 
parallel SWCNTs, and, interestingly, somewhat larger. The dependence of the pair interaction on separation is cx 
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FIG. 2: The van der Waals interaction free energy for two [5,1] SWCNTs (panel a), two [29,0] SWCNTs (panel b) and [5,1] 
and [29,0] SWCNTs (panel c) crossed at the right angle (9 = 7r/2) as a function of their separation. Circles represent the 
nonretarded expression^, while squares represent the fully retarded calculation Eq. 1211 Dashed and dotted lines indicate 
£~ 4 and £~ 5 dependences, respectively. The interaction free energy is calculated only for interaxial separations larger then 
£ — Ri + R2, i.e. in the non-intersecting regime of the cylinders. Close to this value the thin cylinder approximation breaks 
down and is superseded by a different type of calculation, see^ for details. 



for the nonretarded case, but one can see a gradual transition towards the oc £~ 5 dependence for large distances. 
Again, the interaction is strongest for a pair of [29,0] carbon nanotubes. 

The dependence of the van der Waals interaction on the relative orientation angle (9) is of interest, since it can be 
used to calculate the effective van der Waals torque that acts to make the cylinders parallel. The dependence of the 
van der Waals interaction on the angle for I = 4 nm is shown in Fig. [3] Note that the torque vanishes for 6 = ir/2. 
One should also note that for 9 = the van der Waals interaction calculated from Eq. rXU] diverges. This is due to the 
fact the the interaction free energy for parallel cylinders scales with their length and thus diverges for infinitely long 
cylinders. It is clear from Eqs. [T0l and [T2l that the angular dependence of the van der Waals interaction depends on 
the details of the dielectric response of the cylinders. Nevertheless, a quick look at panel a) of Fig. [3]may lead one to 
think that the angular dependencies for all three cases are completely the same up to some multiplicative constant. 
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FIG. 3: Panel a): The van der Waals interaction energy between two two [5,1] SWCNTs (pluses), two [29,0] SWCNTs (x-es) 
and [5,1] and [29,0] SWCNTs (stars) separated by £ = 4 nm as a function of their relative angle (fully retarded calculation). 
Panel b): The same as in panel a), only the energies have been rescaled so that the three cases yield the same value for 6 = 7r/2. 



This is in part due to the divergence introduced by the 1/ sin 9 factor which tends to screen the fine details of the 
interaction. However, a closer look indeed shows that the angular anisotropies are not scalable as is shown in panel 
b) of Fig. [3l where all three interactions have been scaled so as to reproduce the value of [5,1] - [29,0] interaction at 
6 = 7r/2 (this means that the scaling factor for [5,1] - [29,0] interaction is 1). One can see how the curves separate 
as 8 diminishes, and that [5,1] - [5,1] case has the largest angular anisotropy of the van der Waals interaction. The 
effect is relatively small, however, and the l/sinfi 1 factor determines the anisotropy to the largest extent, at least for 
the cases considered here. 
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IV. THE INFLUENCE OF DIELECTRIC RESPONSE OF BATHING MEDIUM ON THE VAN DER 
WAALS INTERACTION: REPULSIVE, ATTRACTIVE, OR BOTH? 



An interesting question is whether the cylinder-cylinder interaction can be repulsive and how that depends on 
the dielectric responses of the cylinders and the bathing medium (this can be achieved only when the cylinders are 
different). The answer is somewhat hidden in Eq. [2T]and much more explicit in Eq. [23l In both of these equations, the 
van der Waals interaction depends on the product Aj mAq y which appears under the sum over Matsubara frequencies. 
There is, however, the additional term that multiplies it [3 + 5(ai+a2) + 19aia2 for parallel cylinders in the nonretarded 
regime], that depends on both the longitudinal and transverse polarizabilities of the cylinders. A closer examination 
of Eq. [6] shows that for isotropic cylinders a is neccessarily positive, so that the aditional term is also positive when 
both cylinders are isotropic, at least in the nonretarded regime. Parameter a is negative only when the dielectric 
response of the medium is between the longitudinal and transverse responses of the cylinder, and even then, the sign 
of the additional term depends on the details of response of other cylinder. It thus seems that the repulsive nature of 
van der Waals interaction can be most easily obtained and comprehended when it is enforced through the longitudinal 
response of the two cylinders. Assuming now that the terms that multiply AjiiAa 11 (two-dimensional integral in the 
retarded case and relatively simple combination of transverse and logitudinal dielectric responses in the nonretarded 
case) are positive, one concludes that the repulsive van der Waals interaction can be obtained by making the product 
A 1 ||A 2i || negative for all imaginary frequencies, i.e. that (see Eq. BJ 



i(iu>) — e m (iu))][e2 ii(iw) — e m (ioj)} < 0, (26) 



(a completely analogous finding has been recently verified experimentally for the repulsive interaction between gold 
sphere and silica plate immersed in bromobenzenc 16 ). This leads us to an interesting conception of "designing" the 
van der Waals interaction between the cylinders by the introduction of an appropriate intervening medium. A simplest 
way which should produce repulsive interaction (at least in the Ax nA2.11 "channel") is to "construct" the medium 
dielectric response as 

£m M = MV 2 2 ' I|V ' . (27) 

Since dielectric responses are difficult to measure experimentally the "designer" response above would have to be 
engineered by the use of ab initio codes to fill this gap. 

The assumption that the above dielectric response of the intervening medium leads to repulsive dispersion interac- 
tions between two SWCNTs can be easily checked numerically and, as shown in Fig. 01 the van der Waals interaction 
between [5,1] and [29,0] SWCNTs is indeed repulsive, at small distances £, but an interesting effect takes place when 
£ w 35 nm. Namely the retarded van der Waals interaction changes sign, becoming attractive when £ > 35 nm and 
with an extremely shallow minimum (~ 10 -10 fcBT/nm) at £ ~ 36 nm. This effect is not seen in the nonretarded van 
der Waals interaction which is repulsive for all distances £. Although we aimed at constructing the medium response 
that shall guarantee the repulsive van der Waals interaction, we obtained that the retardation effects can act so as 
to change the sign of interaction at a certain separation distance. This effect is quite similar to what has been found 
in Ref . [jja and it is of interest to see whether it can be " boosted" and brought to smaller separation distances by a 
carefully guided choice of the medium in between the cylinders. 

Retardation acts so as to screen the contribution of higher Matsubara frequencies to the total value of the van der 
Waals interaction. This effect becomes more important as the separation distances increase 7 . It thus seems possible 
to design a dielectric response of the medium so that the retardation screens the repulsive large n contributions to the 
summation in Eq. 1211 switching from repulsive to attractive behavior for some £. One can even imagine a medium 
in which the van der Waals interaction between the cylinders would be attractive for small £'s and repulsive for large 
ts. This shall be illustrated by several examples. 

Figure [5] shows the longitudinal dielectric responses of [5,1] and [29,0] SWCNTs. The two model medium re- 
sponses are indicated by dashed and dotted lines, respectively. The two responses are modeled as e m (n) = 
1 + 2.215 exp(— 0.015rt) (dashed line, model 1) and e m (n) — 1 + 2.305 exp(— 0.015n) (dotted line, model 2), where n is 
the Matsubara frequency index. The two model medium responses are quite similar. Examination of the nonretarded 
variant of the van der Waals interaction (Eq. |2"3"|) , suggests that there exists a possibility of observing interesting effects 
in the longitudinal "channel" of dielectric response. Namely, the product A^y A 2) || is positive when n < 10, negative 
when 10 < n < 200, and again positive when n > 200. Thus, the contributions to the total (summed) van der Waals 
interaction are both positive and negative, depending on the value of n. As retardation screens the contribution from 
larger values of n at large separation distances, one may a priori expect to see a change in character of van der Waals 
interaction (repulsive vs. attractive) depending on the separation distance. This is indeed confirmed by numerical 
results shown in Fig. [5J However, the overall behavior of the van der Waals interaction, being attractive for model 1 



11 



10 L 



£ 10- 
c 



CO 



CD 



10 c 



10 



12 



retarded 



[5,1] -[29,0] 



nonretarded 



o 

□ o 
□ o 
□ o 



o 

□ o 

□ o 



□ 



□ □ 



□ □ o 

□ o 

D □ ° 

repulsive attractive 



10 

f[nm] 



100 



FIG. 4: The absolute value of the retarded (squares) and nonretarded (circles) van der Waals interaction energy between 
parallel [5,1] and [29,0] SWCNTs as a function of their separation I. The intervening medium is described as Eq. [27] The 
retarded van der Waals interaction is repulsive when i < 35 nm and attractive when £ > 35 nm. 



medium, and repulsive at short distances for model 2 medium, is even more striking. This change can occur by quite 
a minor alteration of the medium response. Note, however, that the absolute magnitude of [5,1] - [29,0] van der Waals 
interaction is a factor of 200 smaller with respect to the case when the SWCNTs are in vacuum (~ 0.005 fc^T/nm in 
medium vs. ~ 1.1 fcflT/nm in vacuum at I = 1.32 nm). An expected minimum in the van der Waals interaction is 
observed in model 2 of the medium at I « 17 nm, but its depth is quite small (~ 10 _8 fcs2V nm )- The origin of this 
effect is exactly the same as in the case studied by Elbaum and Schick - the medium dielectric response (water in their 
case) is such that the product [ej y(iw) — e m (iuj)][e2 ii(iw) — ^(iw)] is both positive and negative, depending on the 
Matsubara frequency, icu n (there is a crossover (or several crossovers) between the medium dielectric response and one 
of the cylinders' responses). In the case studied by Elbaum and Schick, the medium response (water) is quite similar 
to the response of one of the two half-spaces (ice). This leads us to experiment a bit more with the medium response, 
making it more similar to the response of one of the cylinders. The results of one such experiment are shown in Fig. [7] 
The minimum in now located at quite small separations (3.5 nm), and its depth is larger (5.3 10 -6 fc^r/nm). We also 
observe a maximum (~ 2 10~ 4 fc^T/nm at ~ 1.5 nm) in a narrow repulsive region of interaction separating the two 
attractive interaction regions (I < 1.4 nm and I > 2.7 nm). This interesting effect was not observed in the previous 
models. We should also note here that the appearance of the minimum in the interaction is a quite delicate effect 
that needs precise tuning of the parameters of the medium response. From our numerical experiments we found that, 
while it is in general easy to obtain purely repulsive or purely attractive interactions, the minimum in the interaction 
appears as a quite special effect, and only in a narrow region of parameters describing the response of the intervening 
medium. 



V. DISCUSSION 



The derivation of van der Waals - dispersion interaction between two anisotropic dielectric cylinders presented here, 
based on the Pitaevskii ansatz, is valid when the interaction energy either scales linearly with the length, as in the 
case of parallel cylinders, or does not scale with the length at all, as in the case of skewed cylinders. All the cases 
where the finite cylinder length effects come into play arc not covered by this derivation. In addition the cylinder 
radius has to be the smallest length involved in the problem. Another drawback of this method is that the cases of 
infinitely large dielectric response as in the metallic SWCNTs that show a Drude-fike response at zero frequency can 
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FIG. 5: Longitudinal dielectric responses of [5,1] (pluses) and [29,0] (x-es) SWCNTs as a function of the Matsubara frequency 
index. Two model medium responses are indicated by dashed and dotted lines (see text). 
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not be dealt with within the framework proposed here. The zero frequency term in the Matsubara summation would 
have to be treated differently^ for those cases. 

Despite these drawbacks the present approach yields the final result for the calculation of van der Waals dispersion 
interactions between infinitely long semiconducting anisotropic SWCNTs. We derived all the appropriate limits and 
showed how the detais and peculiarities of the dielectric response of various SWCNTs effect this interaction and can 
in general lead also to non-monotonic van der Waals dispersion interactions of a type analogous to those found by 
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Elbaum and Shick in the case of an aqueous layer on the surface of icei^. 

We analyzed in detail also the effect of retardation at sufficiently large inter cylinder spacings. The general conclusion 
is that these effects kick in at separation of about 100 nm for parallel cylinders and at about somewhat smaller 
separation of 50 nm for skewed cylinders. In its retarded form, that decays one distance power faster than the non- 
retarded form, the dispersion interactions between cylinders are in general small, only a fraction of fc^T per nm length 
of the cylinder. Nevertheless these effects are real and can be significant for sufficiently long carbon nanotubes. 

Since the dispersion interactions depend in a complicated way on the dielectric properties of the interacting 
anisotropic cylinders there exists a possibility that between two different types of SWCNTs, with sufficiently different 
dielectric spectra,, the interaction would become repulsive, or even show a non-monotonic separation dependence. 
Such variation of van der Waals - dispersion interaction has been observed in other contexts-^. These delicate effects 
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in themselves do not depend on the fact that the interacting cylinders are anisotropic. In fact, they persist even when 
the cylinders are made of isotropic dielectric materials. 
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